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Convergence  of  a  PI  Coordination  Protocol 
in  Networks  with  Switching  Topology  and  Quantized  Measurements 

Enric  Xargay,  Ronald  Choe,  Naira  Hovakimyan,  and  Isaac  Kaminer 


Abstract —  This  paper  analyzes  the  convergence  properties  of 
a  distributed  proportional-integral  protocol  for  coordination  of 
a  network  of  agents  with  multiple  leaders,  dynamic  information 
flow,  and  quantized  measurements.  We  show  that  the  integral 
term  of  the  protocol  allows  the  follower  agents  to  ‘learn’  the 
reference  rate,  rather  than  have  it  available  a  priori,  and  also 
provides  disturbance  rejection  capabilities. 

I.  Introduction 

Worldwide,  there  has  been  growing  interest  in  the  use  of 
autonomous  vehicles  to  execute  complex  missions  without 
constant  supervision  of  human  operators.  A  key  enabling 
element  for  the  execution  of  such  missions  is  the  availability 
of  advanced  strategies  for  cooperative  motion  control  of 
autonomous  vehicles.  In  [1],  for  example,  the  authors  address 
the  development  of  robust  strategies  for  cooperative  missions 
in  which  a  fleet  of  UAVs  is  required  to  follow  collision-free 
paths  and  arrive  at  their  respective  final  destinations  at  the 
same  time.  The  distributed  protocol  used  for  group  coordi¬ 
nation,  which  was  first  introduced  in  [2],  has  a  proportional- 
integral  (PI)  structure  in  which  each  agent  is  only  required 
to  exchange  its  coordination  state  with  its  neighbors,  and  the 
constant  reference  rate  is  only  available  to  a  single  leader. 
The  integral  term  in  the  consensus  algorithm  allows  the 
follower  UAVs  to  ‘learn’  the  reference  rate  from  the  leader. 

A  generalization  of  this  PI  protocol  was  proposed  in  [3], 
where  the  authors  developed  an  adaptive  algorithm  to  recon¬ 
struct  a  time-varying  reference  velocity  that  is  available  only 
to  a  single  leader.  The  paper  used  a  passivity  framework  to 
show  that  a  network  of  nonlinear  agents  with  fixed  connected 
topology  asymptotically  achieves  coordination.  The  work 
in  [4]  also  used  a  (discrete-time)  PI  protocol  to  synchronize 
networks  of  clocks  with  fixed  connected  information  flow.  In 
this  application,  the  integral  part  of  the  controller  was  critical 
to  eliminate  the  different  initial  clock  offsets. 

This  paper  modifies  the  PI  protocol  in  [1],  [2]  to  include 
multiple  leaders,  and  analyzes  the  convergence  properties  of 
the  protocol  for  coordination  of  a  network  of  agents  with 
dynamic  information  flow  and  quantized  measurements,  a 
topic  that  has  received  increased  attention  in  recent  years  [5]— 
[10].  On  one  hand,  the  use  of  multiple  leaders  in  the 
protocol  improves  robustness  to  a  single-point  failure.  On 
the  other  hand,  the  use  of  finite-rate  communication  links 
and/or  coarse  sensors  motivates  the  interest  in  quantized 
consensus  problems.  The  main  contribution  of  this  paper  is 
twofold.  First,  we  present  lower  bounds  on  the  convergence 
rate  of  the  collective  dynamics  as  a  function  of  the  number 
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of  leaders  and  the  quality  of  service  (QoS)  of  the  network, 
which  in  the  context  of  this  work  represents  a  measure  of  the 
level  of  connectivity  of  the  dynamic  graph  that  captures  the 
underlying  network  topology.  And  second,  we  analyze  the 
existence  of  equilibria  as  well  as  the  convergence  properties 
of  the  collective  dynamics  under  quantized  feedback. 

The  paper  is  organized  as  follows.  Section  II  describes 
the  problem  formulation.  Section  III  presents  the  PI  protocol 
adopted  in  this  paper  and  analyzes  its  convergence  properties. 
In  Section  IV,  we  study  the  collective  dynamics  under 
quantization.  Simulation  results  are  presented  in  Section  V, 
while  Section  VI  summarizes  concluding  remarks. 

II.  Problem  Formulation 

Consider  a  network  of  n  integrator-agents 

Xiit)  =  Ui{t)+di,  xfl 0)  =  xi0,  i6l„:=  {1,..., n},  (1) 

with  dynamic  information  flow  Go(t )  :=  (Vo,£o(i)).  In  the 
above  formulation,  xflt)  £  R  is  the  coordination  state  of  the 
?'th  agent,  uflt)  £  R  is  its  control  input,  and  d,  £  R  is  an 
unknown  constant  disturbance. 

The  control  objective  is  to  design  a  distributed  protocol 
that  solves  the  following  coordination  problem : 

xflt)  -  Xj  ( t )  *^5?  o  ,  V  i,  j  £  ln  ,  (2a) 

xflV^p,  Vieln,  (2b) 

where  p  is  the  desired  (constant)  reference  rate. 

The  network  and  the  communications  between  agents 
satisfy  the  following  assumptions: 

Assumption  1:  The  ?'th  agent  can  only  exchange  informa¬ 
tion  with  a  set  of  neighboring  agents,  denoted  by  A/)°(f). 

Assumption  2:  Communications  between  two  agents  are 
bidirectional  ( Go{t )  is  undirected)  and  the  information  is 
transmitted  continuously  with  no  delays. 

Assumption  3:  The  connectivity  of  Go  (t)  at  time  t  satisfies 
the  persistency  of  excitation  (PE)-like  condition 

t+T 

QnL0{T)QndT  >  pin-!  ,  Vf>0,  (3) 

where  Lo(t)  £  R”x"  is  the  piecewise-constant  Laplacian  of 
the  graph  Go(t),  and  Qn  is  any  (n  —  1)  x  n  matrix  satisfying 
Qnln  =  0  and  QnQn  =  In-i,  with  ln  being  the  vector 
in  R”  whose  components  are  all  1.  The  parameters  T  >  0 
and  p  £  [0, 1]  characterize  the  QoS  of  the  communications 
network,  which  in  the  context  of  this  paper  represents  a  mea¬ 
sure  of  the  level  of  connectivity  of  the  dynamic  graph  Go(t). 

Remark  1:  Condition  (3)  requires  only  the  graph  Go(t)  to 
be  connected  in  an  integral  sense,  not  pointwise  in  time.  In 
fact,  the  graph  may  be  disconnected  during  some  interval  of 


time  or  may  even  fail  to  be  connected  at  all  times.  Similar  The  lemma  below  shows  that  the  connectivity  of  the 
type  of  conditions  can  be  found  in  [11]  and  [12].  graph  Q(t)  satisfies  a  PE-like  condition  similar  to  (3). 


III.  Distributed  Consensus  Protocol 


A.  Addition  of  Virtual  Agents 

The  consensus  protocol  adopted  in  this  paper  introduces 
ni.  virtual  agents  (1  <  nt  <  n )  in  the  network,  associated 
with  rif  agents.  These  virtual  agents  are  implemented  in  nt 
distinct  agents  and  have  the  following  dynamics: 

±u(t)  =uu(t),  xu{ti)  =  xH0, 

where  the  virtual  control  laws  unit),  i  G  It,  are  yet  to  be 
defined.  Without  loss  of  generality,  we  assume  that  these 
virtual  agents  are  implemented  in  agents  1  to  nt,  that  is, 
the  ith  virtual  agent  is  implemented  in  the  /th  agent.  In 
the  context  of  this  paper,  these  nt  agents  are  referred  to 
as  leaders,  while  the  remaining  agents  are  followers. 

To  limit  the  amount  of  information  transmitted  over  the 
network,  each  leader  is  only  allowed  to  exchange  the  state 
of  its  virtual  agent  with  its  neighbors;  in  other  words,  the 
?'th  leader  can  only  transmit  the  state  xu(f),  rather  than 
transmitting  both  Xu (t)  and  Xi(t).  Finally,  we  note  that 
the  agent  and  the  virtual  agent  of  a  leader  can  exchange 
information  uninterruptedly,  as  these  two  agents  do  not 
communicate  over  the  network.  Figure  1  presents  an  example 
illustrating  the  addition  of  two  virtual  agents  in  a  network 
of  three  agents. 


0 


0" 


0 


Lemma  1:  Consider  a  network  with  n  agents  and  rif  vir¬ 
tual  agents,  added  according  to  the  description  above.  If  the 
connectivity  of  the  original  network  satisfies  Assumption  3, 
then  the  connectivity  of  the  extended  network  verifies 

1  1  rt+T 

QNL(r)QTNdT  >  fj,ne Ijv— l,  V  t  >  0, 

where  Qn  is  any  (N  —  1)  x  N  matrix  such  that 
QjvIjv  =  0  and  QnQJj  =  Ijv-ii  while  the  constant  pni 
characterizes  the  QoS  of  the  extended  network.  The  para- 
mater  jint  can  be  determined  recursively  from  the  relation 

/A  =  F(n  +  i  -  1,  A**— i)  ,  i  =  l,...,ne, 
together  with  the  initial  condition  po  =  ft,  and 

LPH  „  -i  _  (fci/j+2)— (fexM+2)2-4x^(fc-|-l) 
t  2(fc+l)  ’ 

Proof.  The  proof  is  omitted  due  to  space  limitations.  □ 


Remark  2:  The  concept  of  virtual  leader  is  quite  common 
in  the  consensus  literature,  where  it  is  understood  as  an 
additional  agent  running  in  open  loop  and  providing  an 
external  consensus  reference  state  to  a  subgroup  of  fol¬ 
lower  agents;  see  [13]— [16]  and  references  therein.  In  our 
framework,  the  virtual  agents  (or  virtual  leaders )  play  a 
different  role  in  the  network  protocol  and  their  dynamics 
are  affected  by  other  agents  (virtual  and  non-virtual).  From 
a  functional  perspective,  the  virtual  agents  in  our  framework 
are  equivalent  to  the  “subgroup  leaders”  in  the  work  reported 
in  [15].  The  key  reason  for  implementing  these  virtual  agents 
as  part  of  our  consensus  protocol  is  to  provide  a  set  of  agents 
with  disturbance-free  dynamics,  which  -as  will  become  clear 
later-  allow  to  effectively  solve  the  coordination  problem  (2). 


(a)  Original  network. 


(b)  Network  with  two  virtual  agents. 


Fig.  1.  Addition  of  =  2  virtual  agents  in  a  network  of  n  =  3  agents. 


The  inclusion  of  these  rif  virtual  agents  results  in  a  new  ex¬ 
tended  network  of  N  :=  n  +  ng  agents  with  a  new  dynamic 
topology  (/(f).  According  to  the  description  above,  this  new 
topology  is  characterized  by  the  following  neighboring  sets: 

Mi  :=  {£i}  ,  i  G  If  , 

Mi  (f )  :=  (A/)0  (f)  \  If)  U  Li  (f) ,  i$lt, 

Mu  (f )  :=  (Mf  (f )  \  h )  U  Ci  (f )  U  {* }  ,  %€Xt, 

where  the  vertex  set  £,;(f)  is  defined  as 

:  j  G  nl<)}  . 


The  Laplacian  L(t)  of  the  new  extended  graph  with  vertex 
set  V  :=  {£1, . . . , Int,  1 ,n}  is  given  by 

L(f)  =  P/  [g  L°w  }Pe  +  Lv  G  M™  , 
where  P#  is  the  (0, 1  ^-permutation  matrix 


Pr  ■■= 


o  i„ 

In*  0 

0  0 


while  Lv  is  defined  as 


-In*  0 
In*  0 

0  0 


p NxN  . 


r,NxN 


B.  Proportional-Integral  Protocol 


To  solve  the  consensus  problem  (2),  we  adopt  the  protocol 


UU  =  kp  (xi  -  xti)  +  Pi  iG  (4) 

u*  =  kP  (XI  ~  0)  +X*’  ieXr 1,  (5) 

Xi  =  ki  (0  “  0)  >  Xi(0)  =  Xm ,  i  G  In,  (6) 

where  kp  >  0  and  kj  >  0  are  coordination  gains.  This 
protocol  has  a  PI  structure  in  which  each  agent  is  only 
required  to  exchange  its  coordination  state  x,(t)  with  its 
neighbors,  and  the  reference  rate  p  is  only  available  to  the 
nt  leaders.  We  also  note  that  the  virtual  agents  adjust  their 
dynamics  according  to  information  exchanged  with  their 
neighboring  agents. 

The  protocol  (4)-(6)  can  be  rewritten  in  compact  form  as 


u(t) 

x(t) 


— kpL(t)x(t )  + 


Pl'n.£ 
x(t)  \  ’ 


-k!CT L(t)x{t) , 


(7) 

X(0)  =  Xo  , 


where  u(t),  x(t),  and  x(f)  are  defined  as 

u(t)  :=  [un{t),  •  ■  ■ ,  utn(t),  iti(f), . . . ,  un(t)\T  G  RN  , 
X (f)  :=  [xtlif),  .  .  .  ,X£n{t),Xi{t),  . .  .,xn(t)\T  G  Rn  , 
x{t)  ■■=  [xi(f),  •  ■  • ,  Xn(f)]T  g  R" ; 
and  CT  :=  [  0  I„  ]  G  RnxN. 


C.  Collective  Dynamics  and  Convergence  Analysis 


Protocol  (7)  leads  to  the  closed-loop  collective  dynamics 


x{t)  =  — kpL(t)x(t ) 


Pt-ng 

x(*)+dj  ’ 


tc(0)  =  Xq  , 


X(t)  =  -k!Cr L(t)x(t) ,  x(0)  =  Xo, 

where  d  :=  [di, . . . ,  dn]T  £  Kn  is  the  disturbance  vector. 
Note  that  the  solutions  (in  the  sense  of  Caratheodory  [17]) 
of  the  collective  dynamics  above  exist  and  are  unique,  since 
the  Laplacian  L(t)  is  piecewise  constant  in  t. 

To  analyze  the  convergence  properties  of  the  algorithm  (7), 
we  reformulate  the  consensus  problem  (2)  into  a  stabilization 
problem.  To  this  end,  we  define  the  projection  matrix  II jv  as 


UN:=lN-^lk  £  RNxN  , 
and  note  that  the  following  equalities  hold: 

ii/v  =  n  N  =  iij^ ,  Q'nQn  =  n 

L(t)TA.iv  —  II^Y L(t)  =  L(f) . 


jv  , 


Moreover,  we  have  that  the  spectrum  of  the  matrix 
L(t)  :=  QNL(t)QTN  £  R(v-i)x(jv-i) 

is  equal  to  the  spectrum  of  the  extended  Laplacian  L(t) 
without  the  eigenvalue  Ai  =  0  corresponding  to  the  eigen¬ 
vector  ljv-  Finally,  we  define  the  consensus  error  state 
C(t)  :=  [CiTW,C2W  as 

Ci  it)\=QNx(t)  £  R"-1 , 

C2{t)  :=  x(t)  -  pin  +  d  £R". 


Note  that,  by  definition,  Ci  (f)  =  C2  (i)  =  0  is  equivalent  to 
x{t)  £  spanjljv}  and  x(t)  =  pljv- 

With  the  above  notation,  the  closed-loop  collective  dynam¬ 
ics  can  be  reformulated  as  (see  Appendix) 


c(t)  =  Ac(t)c(t),  C(o)  =  Co , 

where  A^(t)  £  R(Ar+«-1)x(JV+ra-1)  js  given  by 


Ac(t) 


— kpL(t )  Q n  C 
-kjC^Q^Lp)  0 


(8) 


Next  we  show  that,  if  the  connectivity  of  the  graph  Go(t) 
verifies  the  PE-like  condition  (3),  then  protocol  (7)  solves  the 
consensus  problem  (2).  The  next  theorem  proves  this  result. 


Theorem  1:  Consider  the  collective  dynamics  (8)  and  sup¬ 
pose  that  Go(t)  verifies  the  PE-like  condition  (3)  for  some 
parameters  p  and  T.  Then,  for  any  kp  >  2,  there  exist 
coordination  gains  kp  and  kq  such  that  the  inequality 

||C(t)ll<«cllC(0)||  e~Xet 

holds  for  some  positive  constant  a £  G  (0,oo),  and  with 


A/i  ^  A/>  .  — 


kpN  fin 


(!  +  **£) 


JV  \-l 


(1  +  fcpJVT)2  '  --Pnt‘ 

Also,  the  coordination  states  and  their  rates  of  change  satisfy 

limt-j-oo  | Xi(t)  -  Xj(t) |  =  0 ,  i,j  £  Tn  , 

limt-j-oo  Xi(t)  =  p,  i  G  ln- 


Proof.  The  proof,  which  is  omitted  here  due  to  space 
limitations,  is  similar  to  the  proof  of  Lemma  3  in  [1].  □ 


Remark  3:  Theorem  1  above  indicates  that  the  QoS  of 
the  network  (characterized  by  T  and  p)  limits  the  achievable 
(guaranteed)  rate  of  convergence  of  the  closed-loop  collective 
dynamics.  According  to  the  theorem,  for  a  given  QoS  of  the 


network,  the  maximum  (guaranteed)  rate  of  convergence  A* 
is  achieved  by  setting  kp  =  which  results  in 

Ac*:=^  (f +kf)%y1 . 

We  also  note  that,  as  T  goes  to  zero  (graph  connected 
pointwise  in  time),  the  convergence  rate  can  be  set  arbitrarily 
fast  by  increasing  the  coordination  gains  kp  and  kj. 

IV.  Convergence  under  Quantization 
In  this  section  we  analyze  the  stability  and  performance 
characteristics  of  the  distributed  PI  protocol  presented  in 
the  previous  section  when  the  agents  exchange  quantized 
measurements.  For  the  sake  of  simplicity,  in  this  paper  we 
consider  only  uniform  quantizers  with  step  size  A. 


A.  Protocol  and  Collective  Dynamics 

When  only  quantized  information  from  the  other  agents  is 
available,  the  PI  protocol  introduced  in  (7)  becomes 

-kP  (. D(t)x  -  A{t)  q(a:))  +  [pl™‘  ]  , 


U  =  —h 


X  =  -kiC] 


(9) 


^D(t)x  -  A(t)  q(a:)J  ,  x(0)  =  Xo, 
where  q (x(t))  £  ZWA  is  the  quantized  coordination  state 
q (x{t))  :=  [qA(xa(f)), . . . ,  qA (xn(t))]r  , 
with  qA(-)  :  R  —>  ZA  being  defined  as 

qA(0  :=  sgn(£)A  |f  +  |J  ,  ££R. 

The  time- varying  matrices  D{t)  and  A(t)  are  defined  as 
D(t)  :=  D(t)  +  D(, ,  A{t)  :=  A(t)  +  D£  , 

where  D(t)  and  A(t)  are  respectively  the  degree  and  adja¬ 
cency  matrices  of  L(t),  while  D(  is  given  by 


Dr  := 


0 

-In 


nJVx  JV 


Note  that  only  the  information  exchanged  over  the  network 
is  subject  to  quantization;  in  fact,  each  agent  has  access  to 
its  own  unquantized  state,  and  leaders  also  have  access  to 
the  unquantized  state  of  its  virtual  agent  (and  viceversa). 

Then,  noting  that  L(t)  =  D(t)  —  A(t),  the  collective  dy¬ 
namics  can  be  written  as 


x  =  — kpL{t)x  + 


Pt-Tl£ 

X+d 


+  kPA(t)ex  ,  tc(0)  =  x0  , 


X  =  -kICTL(t)x  +  kjCTA(t)ex  ,  x(0)  =  Xo  , 


where  ex(t)  :=  q(tc(f))  —  x(t)  is  the  quantization  error 
vector.  In  terms  of  the  consensus  error  state  C(i),  the 
collective  dynamics  can  be  expressed  as 

CO)  =  A;(£K(f)  +  Bc(t)ex(t),  C(0)  =  Co,  (io) 


where  A^{t)  was  introduced  in  (8)  and  B^(t)  is  given  by 


BC(t)  := 


fepQjv  A(t) 
fc,CTA(t) 


Note  that,  in  this  case,  the  right-hand  side  of  the  collective 
dynamics  is  discontinuous  not  only  due  to  the  time-varying 
topology,  but  also  due  to  the  presence  of  quantized  states. 
As  proven  in  [8],  Caratheodory  solutions  might  not  exist 
for  quantized  consensus  problems,  implying  that  a  weaker 
concept  of  solution  has  to  be  considered.  Similar  to  [8],  we 
will  consider  solutions  in  the  sense  of  Krasovsky  [17]. 


To  show  that  Krasovsky  solutions  to  (10)  exist  (at  least) 
locally,  we  note  that,  during  continuous  evolution  of  the 
system  between  “quantization  jumps”,  the  network  dynam¬ 
ics  (10)  are  linear,  with  the  quantized  state  q (x(t))  acting 
as  a  bounded  exogenous  input.  This  implies  that  the  solu¬ 
tions  x(t)  are  locally  bounded  (no  finite  escape  time  occurs). 
Then,  local  existence  of  Krasovsky  solutions  is  guaranteed 
by  the  fact  that  the  right-hand  side  of  (10)  is  measurable  and 
locally  bounded  [17].  At  this  point,  we  cannot  claim  that 
Krasovsky  solutions  to  (10)  are  complete;  for  this,  we  will 
need  to  prove  that  solutions  are  bounded  (see  Theorem  2). 

B.  ( Krasovsky )  Equilibria 


error  state  degrades  gracefully  with  the  value  of  the  quantizer 
step  size.  The  next  theorem  summarizes  this  result. 

Theorem  2:  Consider  the  closed-loop  collective  dynam¬ 
ics  (10)  and  suppose  that  the  topology  Go(t)  verifies  the  PE- 
like  condition  (3)  for  some  parameters  //  and  T.  Then,  there 
exist  coordination  gains  kp  and  ki  ensuring  that  there  is  a 
finite  time  Tf,  >  0  such  that  the  bounds 

I Xi{t)  -  Xj{t) |  <  avA  , 

| Xi{t)  —  p\<  apA  , 

hold  for  all  t>Tt  and  some  constants  av,ap  €  (0,  oo). 
Proof.  The  proof  is  omitted  due  to  space  limitations.  □ 


Before  investigating  the  convergence  properties  of  the 
quantized  collective  dynamics  (10),  in  this  section  we  an¬ 
alyze  the  existence  of  equilibria  for  these  dynamics.  To 
simplify  the  analysis,  we  assume  (only  in  this  section)  that 
the  network  topology  is  static  and  connected.  Under  this  as¬ 
sumption,  one  can  easily  show  that  the  unquantized  collective 
dynamics  (8)  have  one  isolated  equilibrium  point  at  £eq  =  0. 
However,  when  quantized  information  is  exchanged  over  the 
network,  £eq  =  0  is  not  an  equilibrium  point  of  the  collective 
dynamics  anymore  and  other  (undesirable)  equilibria  might 
exist,  depending  on  the  step  size  of  the  quantizers. 

To  show  this,  we  first  notice  that  £(i)  =  0  is  equivalent 
to  x(t)  £  spanjljv}  and  x(t)  =  0  holding  simultaneously. 
Hence,  £eq  :=  [Ci'eq-  Ciq]  is  an  equilibrium  of  (8)  if 


7(f)ljV  £  1C  (—kp  (^Dxeq(t)  —  jiq(leq(l))j  +  £«+d  ) 


0  g  AC  - 


^  kjC  ^JDa?eq(f)  A  q(a?eq 


where  7  (t)  €  R  is  an  arbitrary  signal;  tceq(f)  is  a  continuous 
coordination-state  trajectory  satisfying  (jieq  =  Qjv*eq(i); 
while  Xeq  :=  C2eq  —  pi™  +  d.  The  second  inclusion  above 
and  continuity  of  a;eq(f),  along  with  the  fact  that  the  network 
is  assumed  to  be  static  and  connected,  preclude  the  existence 
of  equilibria  involving  time-varying  coordination-state  trajec¬ 
tories,  i.e.  7 (f)  =  0  (or  equivalently  £ceq(f )  =  0).  Then,  the 
set  of  (Krasovsky)  equilibria  of  (8)  can  be  defined  as 


E  : —  <  (iEeq,  Xeq)  6 


r,N 


[  j 

Xeq+dJ 


(ID 


Q  j£  (  — fcp(Das0q-Aq(a!Gq))+^ 

V  —  kiCT  (T)a:eq  —  A  q(a:eq)) 

Next,  we  show  that,  under  sufficiently  fine  quantization, 
the  set  E  is  empty. 


Lemma  2:  Consider  the  quantized  collective  dynam¬ 
ics  (10),  and  assume  the  network  topology  is  static  and 
connected.  If  the  step  size  of  the  quantizers  satisfies 

A  <-  2t i,  |p| 
n(n— 1)  kp  5 


(12) 


then  the  set  of  equilibria  E  is  empty. 

Proof.  The  proof  is  omitted  due  to  space  limitations.  □ 


C.  Convergence  Analysis 

Next  we  show  that,  if  the  connectivity  of  Go  (t)  verifies  the 
PE-like  condition  (3),  then  protocol  (9)  solves  the  consensus 
problem  (2)  in  a  practical  sense.  Moreover,  the  consensus 


V.  Simulation  Results 

We  now  present  simulation  results  illustrating  the  theo¬ 
retical  findings  of  the  paper.  To  this  end,  we  consider  a 
network  of  5  agents  with  dynamics  (1).  At  a  given  time  t, 
the  information  flow  is  characterized  by  one  of  the  graphs 
in  Figure  2;  note  that  all  four  graphs  are  not  connected.  The 
control  objective  is  to  design  a  distributed  PI  protocol  that 
solves  the  consensus  problem  (2)  with  p  =  1  (in  a  practical 
sense).  In  all  of  the  simulations,  the  initial  coordination-state 
vector  xo  and  the  disturbance  vector  d  are  given  by 

XQ  =  [-1,  2,  4,  -4,  3  ]T  ,  d=[  o,  5,  -3,  4,  l]T  . 

To  solve  the  consensus  problem,  we  add  2  virtual  agents 
to  the  network,  and  implement  the  (quantized)  protocol  (9) 
with  PI  gains  kp  =  0.60  and  fcj  =  0.15,  and  initial  integrator 
state  xo  =  0.  Figure  3  presents  the  computed  evolution 
of  the  closed-loop  collective  dynamics  with  quantizer  step 
size  A  =  0.3  (note  that  this  step  size  verifies  inequality  (12)). 
The  figure  shows  the  time  evolution  of  the  coordination 
states,  their  time-derivative,  the  integrator  states,  and  the 
2- norms  of  the  consensus  error  states  £i(i)  and  £2  (i). 
Additionally,  Figure  4  shows  an  estimate  of  the  QoS  of  both 
the  original  network  and  the  extended  network,  computed  as 

p(t)  ■=  Amin  (A  Y  ILt  QnL0(T)Q^dT^j  ,  t>T, 

Pm(t)  ■=  Amin  y  //_T  ,  t>T, 

with  T  =  0.25  sec.  The  results  demonstrate  that  the  PI  dis¬ 
tributed  protocol  allows  the  followers  to  ‘learn’  the  reference 
rate  command  p  and  reach  agreement  with  the  leaders, 
while  effectively  compensating  for  the  (constant)  distur¬ 
bances  present  in  the  network.  To  illustrate  the  effect  of 


(c)  Topology  3  (d)  Topology  4 

Fig.  2.  Network  topologies. 
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(c)  Rate  of  the  coordination  states,  x(t)  (d)  Integrator  states,  xW 

Fig.  3.  Quantized  closed-loop  collective  dynamics  with  fine  quantization  (A  =  0.3). 


(b)  Extended  network 

Fig.  4.  Quality  of  service  (T  =  0.25  sec). 


the  QoS  on  the  convergence  rate  of  the  collective  dynamics, 
Figure  5  presents  convergence  times1  (normalized  to  the  case 
of  2  leaders  and  complete  graph  as  network  topology)  as  a 
function  of  the  parameter  fi  (with  T  =  0.25  sec).  The  figure 
shows  that  the  speed  of  convergence  decreases  with  (z)  the 
QoS  of  the  network  and  (zz)  the  addition  of  virtual  agents. 
The  latter  can  be  explained  by  the  reduction  of  the  QoS  of  the 
extended  network  as  virtual  agents  are  added  (see  Lemma  1). 
Notice  that  these  results  are  consistent  with  Theorem  1. 
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Fig.  5.  Convergence  time  of  the  collective  dynamics  (normalized  to  the  case 
of  2  leaders  and  complete  graph  as  network  topology;  fnorm  =  28.4  sec). 

Next,  we  use  the  same  simulation  scenario  to  illustrate 
the  importance  of  adding  virtual  agents  to  the  network.  For 
this  purpose,  we  consider  the  network  of  5  agents  above  and 
use  protocol  (4)  to  drive  agents  1  and  2,  while  agents  3,  4, 
and  5  are  driven  with  protocol  (5)-(6).  The  key  difference 
with  the  simulation  in  Figure  3  is  thus  that  protocol  (4) 
is  applied  directly  to  the  (uncertain)  agents,  rather  than  to 

1  Convergence  time  is  defined  here  as  the  time  it  takes  for  the  2-norm  of 
the  consensus  error  state  £(t)  to  converge  to  a  2%-tube  of  its  initial  value. 


the  corresponding  (disturbance-free)  virtual  agents,  which  are 
not  implemented  in  this  case.  Figure  6  shows  the  evolution  of 
the  collective  dynamics  under  the  same  information  flow  (see 
Figure  4)  and  with  quantizer  step  size  A  =  0.3.  As  can  be 
seen,  the  agents  do  not  reach  agreement  and  the  coordination 
states  do  not  evolve  at  the  desired  reference  rate  p  =  1. 

The  same  network  of  5  agents  and  information  flow  are 
now  used  to  verify  that  the  multi-leader  PI  protocol  (with 
virtual  agents)  is  robust  to  the  loss  of  a  leader.  In  this  case, 
we  simulate  the  sudden  loss  of  agent  2,  which  is  one  of 
the  leaders.  Figure  7  presents  the  response  of  the  collective 
dynamics,  which  shows  that,  despite  the  loss  of  one  of  the 
leaders,  the  PI  protocol  is  still  able  to  solve  the  coordination 
problem  (2).  We  note  that,  for  the  agents  to  reach  the  desired 
agreement,  it  is  required  that  the  resulting  information  flow 
still  satisfy  a  PE-like  condition  similar  to  (3). 

Finally,  the  same  scenario  is  used  to  illustrate  the  ex¬ 
istence  of  undesirable  attractors  in  the  presence  of  coarse 
quantization.  For  this  purpose,  we  change  the  quantizer  step 
size  to  A  =  3.  The  computed  response  of  the  closed-loop 
collective  dynamics  is  shown  in  Figure  8.  In  this  case,  the 
agents  do  not  reach  the  desired  agreement  and,  in  fact, 
the  solution  converges  to  a  neighborhood2  of  one  of  the 
(Krasovsky)  equilibrium  points  characterized  by  (11). 

VI.  Conclusions 

In  this  paper  we  analyzed  the  convergence  properties  of 
a  distributed  PI  protocol  to  coordinate  a  network  of  agents 
subject  to  constant  disturbances.  We  addressed  the  situation 
where  each  agent  transmits  only  its  coordination  state  to  only 

^Notice  that  the  results  in  Section  IV-B  are  derived  for  network  topologies 
that  are  both  static  and  connected;  instead,  the  simulations  presented  here 
consider  a  time-varying  information  flow. 
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(b)  Consensus  error  states,  HCiMII 
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Fig.  6.  Collective  closed-loop  dynamics  with  no 
virtual  agents  (A  =  0.3). 


Fig.  7.  Robustness  of  the  protocol  to  the  loss 
of  a  leader  (A  =  0.3). 


Fig.  8.  Collective  closed-loop  dynamics  with 
coarse  quantization  (A  =  3). 


a  subset  of  the  other  agents,  as  determined  by  the  network 
topology.  Furthermore,  we  considered  the  case  where  the 
graph  that  captures  the  information  flow  is  not  connected 
during  some  interval  of  time  or  even  fails  to  be  connected  at 
all  times.  We  also  analyzed  the  convergence  properties  of  the 
protocol  when  the  agents  exchange  quantized  measurements. 
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Appendix 

Closed-Loop  Collective  Dynamics 

From  the  definition  of  Ci  (£)  and  C2  (i )  and  the  coordination-state 

dynamics,  it  follows  that 

C'lCO  =  —kpQisrL(t)x(t)  +  pQjvIat  +  Qn  [<a°(t)] 

=  —kpQNL(t)x(t)  +  QivC&it) . 

The  properties  of  the  projection  matrix  IIjv,  along  with  the  fact 

that  QnQn  —  imply  that 

Ci  (t)  =  —  kpQ  nUn  L(t)TlNx(t)  +  QisfCCiit) 

=  —  kpQ nQJjQ n L(t)Q]^Q Nx{t)  +  QisrC&it)  (13) 
=  -kpL{t)Ci(t)  +  QNC^{t). 

Similarly,  it  follows  that 

Ca (t)  =  -kTCT L(t)x{t)  =  -krCT QjfL(t)Ci(t) .  (14) 

Equations  (13)  and  (14)  lead  to  the  dynamics  (8). 


